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TECHNICAL PAPER 


THREE-DIMENSIONAL BAROCLINIC INSTABILITY OF A HADLEY CELL 
FOR SMALL RICHARDSON NUMBER 

I. INTRODUCTION 


Stably-stratified baroclinic flow can be destabilized by several different mechanisms. The most 
extensively studied of these instabilities is usually referred to simply as baroclinic instability. For this 
type of instability the Richardson number, Ri, is much greater than unity, and the quasigeostrophic 
equations [ 1 ] are valid. The perturbations of maximum growth rate have a wavelength in the zonal 
(azimuthal) direction and, in the absence of horizontal shear, the gravest possible structure in the 
meridional direction [2,3]. This instability has been observed in the laboratory and is the cause of mid- 
latitude wave cyclones [4,5]. 

A second type of instability, known as symmetric baroclinic instability, occurs for Ri of order 
unity and requires a more general set of equations than the quasigeostrophic set. In this case, the per- 
turbations of maximum growth rate have meridional structure but no, or only weak, zonal structure. 
The study presented in this paper is concerned with this type of instability. Since this instability occurs 
in only a small range of Ri, there would appear to be few geophysical applications. However, Stone [6] 
conjectured that the banded structure of Jupiter’s atmosphere is caused by this instability, and Bennets 
and Hoskins [7] and Emanuel [8] suggested that rain bands and squall lines may be due to this 
instability. Two of the earliest analyses of symmetric baroclinic instability are those of Solberg [9] and 
Kuo [10], In these papers purely two-dimensional perturbations with no zonal structure were considered. 

Stone [11,12,13] was the first to examine the stability of baroclinic flow for small Ri with 
respect to three-dimensional perturbations. For the basic state Stone took a plane parallel baroclinic flow 
with constant temperature gradients and a constant vertical shear consistent with the thermal wind 
balance. Horizontal shear was excluded. Viscous and thermal diffusion effects were neglected in both the 
basic state and perturbation analyses. In this paper such a basic state will be referred to as the Eady 
basic state. Stone found that conventional baroclinic instability dominates if Ri > 0.95; symmetric 
baroclinic instability dominates if 0.25 < Ri < 0.95; and Kelvin-Helmholtz instability dominates if Ri < 
0.25. For symmetric instability, Stone found that the wavelength of maximum growth rate is zero. 
Stone’s results motivated a search for symmetric baroclinic instability in the laboratory [ 1 4, 1 5 ] , and 
some evidence of the predicted meridional structure was observed. However, Stone’s theoretical model 
was substantially different from the experimental reality, and it cannot be concluded that the existence 
of symmetric baroclinic instability in the laboratory has been established. 

In the theoretical studies on symmetric baroclinic instability which followed Stone’s analyses, 
workers added viscous and thermal effects to make the models more realistic. However, in all but one of 
these studies [16], attention was restricted to two-dimensional perturbations with no zonal structure. 
McIntyre [17] considered an Eady basic state with viscosity and thermal diffusivity included in the 
perturbation analysis. He chose an unbounded model. McIntyre found, in the absence of horizontal shear, 
that for the Prandtl number, o, equal to unity, the critical Richardson number, Ri c , is also equal to 

unity, and for a § 1, Ri c increases above unity. These results are shown in Figure la. Thus, diffusive 

effects further destabilize the flow. McIntyre’s choice of a length scale and unbounded model were such 
that the Ekman number, E, disappears formally from his problem. Emanuel [8] considered the same 



model as McIntyre but for a vertically-bounded flow. Emanutl studied (1) hydrostatic disturbances and 
(2) nonhydrostatic disturbances in a neutrally-stratified fluid. He investigated only neutral modes by 
invoking the principle of exchange of stabilities. For both of the above cases Emanuel found that Ri c 

depends on both a and E. He also established that the wavelength of the most unstable normal mode is 
determined primarily by the depth of the fluid and the slopes of the isentropic surfaces rather than by 
the diffusive properties of the fluid. 

In all of the theoretical studies cited above, the Eady basic state was used. In an experimental 
apparatus whose upper and lower boundaries are stationary, such a basic state is a good approximation 
in the interior of the fluid only; it is not a good approximation near the boundaries especially when E 
is not vanishingly small. Antar and Fowlis [18] (hereinafter referred to as AF) presented a two- 
dimensional theoretical study of symmetric baroclinic instability for a fluid contained between two 
horizontal plates of infinite extent. The analysis used a realistic basic state which was obtained as an 
exact analytical solution of the full set of governing nonlinear equations. This solution included both the 
Ekman and thermal layers adjacent to the boundaries. The stability analysis also included viscous and 
thermal diffusion effects and was performed numerically. Antar and Fowlis found that the instability 
sets in when Ri is close to unity and that Ri is a strong function of t ih a and E. Figure la shows these 
results for RL as a function of a for fixed E and compares their results with those found by McIntyre 

[17]. For fixed a, Antar and Fowlis found that Ri c decreases with increasing E until a critical value of 

E is reached beyond which the flow is stable. These results are shown in Figure 1 b. The analysis was not 
restricted to critical values, and the influence of o and E on the growth rate was also determined. For 
the range of parameter values considered, the most unstable wavelength was about half the depth. The 
nonlinear basic state was also used by Antar and Fowlis [19] in a study of the conventional baroclinic 
instability. More complete discussions of the previous theoretical work dealing with two-dimensional 
studies of symmetric baroclinic instability have been given by Emanuel [8] and AF. 

Recently Busse and Chen [16] extended Stone’s studies by including thermal and viscous diffu- 
sion effects in a three-dimensional analysis of symmetric baroclinic instability. The basic state was an 
Eady basic state in which the horizontal boundaries moved consistent with the thermal wind to eliminate 
Ekman layers on the boundaries. They showed in the limit of small E, that the mode of symmetric 
instability at maximum Ri c exhibits an angle of inclination with the direction of the basic state flow. 

Busse and Chen remarked that owing to this effect the range of Ri for which the instability occurs may 
be increased significantly beyond the limits derived by McIntyre [17]. Busse and Chen obtained then- 
solution through an expansion in the zonal wavenumber; and, thus, their solution is restricted to small 
values of this wavenumber. Also, their analysis was confined to critical values; and, hence, growth rates 
were not discussed. 

The objectives of the present study were to extend the previous work on three-dimensional 
baroclinic instability for Ri ~ 0(1). The work of Stone [11,12,13] is extended to a more realistic 
model involving viscous and thermal dissipation and the work of Busse and Chen [16] to finite growth 
rates and a wider range of the zonal wavenumber. The basic state chosen was that used previously by 
Antar and Fowlis [19] and AF. The stability analysis was performed through numerical integration of 
the full perturbation equations. Critical values and growth rates are presented. 

In Section II the equations and methods of solution for both the basic state and the stability 
analysis are presented. Section III deals with the results of the stability analysis, and the conclusions 
are given in Section IV. 


II. THE EQUATIONS AND THEIR SOLUTION 


A. The Basic State 

A Boussinesq fluid, confined between two horizontal plates which are set a distance, d, apart, 
is considered. The coordinate system used is rectangular Cartesian with axes (x, y, z) corresponding to 
the eastward, northward, and vertical directions, respectively, with the origin set midway between the 
plates. The plates and the fluid are assumed to extend to infinity in the x-direction and to a distance 
L (L » d) in the y-direction. Also, the plates and the fluid are taken to rotate as a whole about the 
vertical axis with a constant angular velocity, £2. To maintain baroclinicity, a temperature distribution 
is imposed on both plates in which the temperature is made to decrease in the y-direction. Also, to 
assure that the vertical stratification, and hence the Richardson number, can be arbitrarily and externally 
fixed, the temperature of both the upper and lower plates are set to differ uniformly by a constant 
amount, AT, for all y. A sketch of the model considered is shown in Figure 2. 

The basic state velocity and temperature fields are governed by the Navier-Stokes, energy, and 
mass conservation equations. These equations, in dimensionless form, for the configuration described 
above, and for a two-dimensional, steady state in a rotating reference frame can be reduced to the 
following: 


-2V = E 3 2 U/3z 2 

(1) 

2U = E 3 2 V/3z 2 - 3p/3y , 

(2) 

V3T/3y = (E /ae) 3 2 T/3z 2 

(3) 

T = 3p/3z , 

(4) 


where Y is the velocity vector (U,V), p the pressure, T the temperature and 
e = agy/Q, 2 
E = p/nd 2 
a = v/k 


are a thermal Rossby number, the Ekman number, and the Prandtl number, respectively. In equations (1) 

through (4), length, time, velocity, and temperature were made dimensionless using d, £2'1, and a gdy/fl, 
and yd, respectively, a is the coefficient of thermal expansion, y is the imposed horizontal boundary 
temperature gradient, and g is gravity. A further assumption used in deriving equations (1) through (4) 




a; 


is that throughout the region of interest in the fluid away from the channel walls, 
velocity component, W, is negligible. 

y < L, the vertical 

The solution to equations (1) through (4) which is consistent with the no-slip and the perfectly 
conducting boundary conditions at the solid boundaries; i.e.. 

G 

II 

< 

II 

o 


> at z = ± 1/2 

(5) 

T = ± (AT/2) - y 


is given by 


U(z) = -f(z)/8 + z/2 

(6a) 

V(z) = ~g(z)/8 

(6b) 

T(y,z) = -y + (AT+oe/4)z - oef(z)/16 

(6c) 


where 

f(z) = [coshR (z+l/2)cosR(z-l/2) - coshR(z-l/2)cosR(z+l/2)]/h(R) , 
g(z) = [sinhR (z+l/2)sinR(z-l/2) - sinhR(z-l/2)sinR(z+l/2)]/h(R) , 
h(R) = sinh 2 (R/2) + sin 2 (R/2) , 

R = E" 1 / 2 . 

An assumption made in obtaining the above solutions is that the dimensionless horizontal temperature 
gradient, dT/dy, is constant and equal to -1 throughout the fluid. The assumptions and limitations of 
this model are discussed more fully in Antar and Fowlis [19]. 

The velocity and temperature profiles given by equation (6) are shown in Figure 3 for repre- 
sentative values of E, e, a, and AT. 
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B. The Perturbation Equations 


To study the stability of the stationary basic state, the dependent variables (temperature, velocity, 
and pressure) are first decomposed into basic state and perturbation components. The equations governing 
the perturbation components are obtained by substituting the variables into the Navier-Stokes and energy 
equations, subtracting the basic state, and linearizing. Since the resulting equations are linear, with 
coefficients depending on z alone, they admit of separable solutions of the form: 


[u',v',w',0',p'] = [u(z), v(z), 0(z), p(z)] exp [i(kx+8y-wt)] 


where n'(u',v',w'), p' and 0' are the perturbation velocity vector, pressure, and temperature, respectively. 

This solution represents a traveling oblique wave with frequency w and wavenumber components 
k and 8 in the zonal and the meridional directions, respectively. After substituting the above solution 
into the governing perturbation equations and eliminating the pressure, the linearized momentum and 
energy equations can be written as 


i[w^(kU+8V)]V 2 w+ie(kD 2 V+BD 2 U)w-2i (kDv-8Du) = -EV 2 V 2 w + m 2 0 (7) 

i[w-e(kU+8V)] [m 2 v-iBDw]-ekw(kDV-8DU)-2ikDw = -EV 2 (m 2 v-i8Dw) (8) 

io[w-e(kU+8V)]0 + eo (v- w DT) = -EV 2 0 (9) 


where 


Dw + iku + i8v = 0 


V 2 = D 2 - m 2 
m 2 = k 2 + 8 2 


and D = d/dz is a differential operator. In the above U, V, and T are the basic state velocity and tem- 
perature fields which are non-simple functions of z. The perturbation velocity and temperature vanish 
on the plates implying that 


w = Dw = v = 0 = 0 at z = ± 1/2 


( 10 ) 


The problem defined by equations (7) through (9) with the boundary conditions (10) define an 
eigenvalue problem of the form 


5 


F (E, e, a, AT, k, 8, w) = 0 


( 11 ) 


where the frequency « is in general complex for the temporal stability problem. The basic state is stable 
or unstable depending on whether Im (co) = coj < 0 or wj > 0, respectively. For marginal stability for 

which ojj = 0, Re (co) = a> r , and either E or e are chosen as the eigenvalues. In most of th" ; • ' 

described in the next section, the growth rate, wj, was used as a function of any of the par: ncters tc 
determine the most unstable mode. 

As indicated in the introduction, previous studies of symmetric baroclinic instability have selected 
the Richardson number as a parameter. However, due to the specific non-dimensionalization chosen for 
this study, Ri does not appear as an explicit parameter in equation (11). Ri is defined as the ratio of 
the vertical stratification to the square of the shear,* i.e., 

/ 3T* \ //au*\ 2 

Rl ’ / (if) <12) 

where the asterisk denotes a dimensional quantity. Since in the present problem this definition implies 
a variable Ri as a function of height, a more convenient definition of Ri is its value at mid-depth in the 
channel. Substituting the values of the temperature and velocity gradients of the basic state [equation 
(6)] at z = 0 into the above definition we obtain 

Ri 0 = 4e' 1 (AT ae/4) + 0(Re' R / 2 ) (13) 


which is a constant defined exclusively by the basic state. Note, that for E < 10* 2 (i.e., R > 10), which 
is the range of E for all the calculations presented in Section III, the last term in equation (13) is neg- 
ligible. Since Ri 0 is a complicated function of the basic state, care is required in selecting values of the 

basic state parameters to achieve desired values of Ri Q . (From this point on the suffix is dropped from 
Ri Q .) Examination of equation (13) reveals that to obtain values of Ri close to unit - , AT must be greater 

or less than zero depending on whether o < 1 or o > 1, respectively. Note that for some values of 
AT < 0, there does exist a large region in the interior of the fluid with stable stratification. This is shown 
in Figure 4 which is based on parameter values for the computations whose results are discussed in 
Section III. This problem of the selection of the basic state parameters to achieve specific values of Ri 
is discussed more fully in AF. 

Note that setting k = 0 in the perturbation equations allows for only axisymmetric waves propa- 
gating in the meridional direction, and setting fi = 0 allows for only waves propagating in the zonal 
direction. 


C. Method of Solution 

The governing differential equations (7) through (9) with boundary conditions [equation (10)] 
form an eigenvalue problem for which a nontrivial solution should be available. Due to the complexity 


6 



a 





of the coefficients of these equations (a consequence of the non-simple basic state), a closed form solu- 
tion cannot be found. The only way to solve this problem for the ranges of all the parameters involved 
is by numerical means. There are basically two direct methods for the numerical s 'ution of differential 
eigenvalue problems. These are the matrix and the shooting methods. Due to the complexity of this 
system (viz, a system of coupled, eight ordinary differential equations), the shooting technique is the 
most straightforward to implement. For small values of the Ekman number, the shooting procedure can 
be successful only if it is used in conjunction with a filtering or orthonormalization scheme (see for 
instance, Conte [20]). Otherwise, the inevitable presence of the rapidly growing solution of the differ- 
ential system can quickly render the linearly independent solution numerical- ' dependent. 

For the solution of the present eigenvalue problem, a computer code was developed with an 
eighth-order, variable step, Runge-Kutta-Fehlberg initial value integrator. For the iteration procedure a 
Newton-Raphson method was used, and the orthonormalization process was implemented at each inte- 
gration step. All of the eigenvalues which are presented in the next section were produced with a rela- 
tive tolerance of 10"^ in the iteration process (i.e., all of the eigenvalues presented are correct up to 
four significant figures). For small enough Ekman numbers, the Ekman and the thermal layers along the 
plates are very thin, and the code was made to take at least 10 steps within each layer to assure adequate 
representation of these layers. Normally, this method requires an initial guess for the eigenvalue, and the 
whole process for obtaining the results is started by first trying several guesses. Once a convergence on an 
eigenvalue was obtained, that eigenvalue was used as the initial estimate for obtaining a neighboring one, 
and the process was repeated. Typically, if an eigenvalue existed, a convergence was achieved in under 
10 iterations. 


III. RESULTS AND DISCUSSION 


In this section the results of a three-dimensional study of baroclinic instability for Ri ~ 0C 
are presented. All of the results were obtained using the basic state defined by equation (6) anr’ 
means of the solution for the complete eigenvalue problem described by equations (7) through 
and represented functionally by expression (11). We begin by presenting plots of the growth rate, u- v 

as a funvtion of the zonal wavenumber, k, for selected values of the meridional wavenumber, £, and for 
a range of values of Ri and a. We chose to vary a since previous work on two-dimensional instability 
(e.g., McIntyre [17]; AF; see Section I) had shown that both Ri c and the growth rates vary with a. 

The actual values of a were chosen to be close to unity and in a range for real fluids. For all the results 

presented, E = 10~ 3 ; this choice was based on the previous work (AF) which showed that for two- 
dimensional instability, Ri c is close to its asymptotic limit for this value. A smaller value would have 

substantially increased the computing time (see Section II). Further, E = 10 3 is a realistic value for 
laboratory experimental studies. The values of AT were selected so that Ri ~ 0(1) (see Section II). 

Figures 5 through 8 show cuj versus k for four values of Ri, o = 2 and AT = -1.5. Figure 5 for 

Ri = 0.93 shows that pure symmetric waves (k = 0) do grow and that their growth rates increase with 
increasing £ reaching a maximum at about £ = 8 and then decrease for a further increase in 8. However, 
the overall maximum growth rate is associated with the conventional baroclinic instability mechanism 
(C = 0). The mode (k = 0.88, 6 = 0) has the maximum growth rate. Figure 6a shows that for Ri = 0.80, 
the growth rates for the symmetric modes have increased substantially and are comparable to the growth 
rates of the purely zonal waves (£ = 0). The maximum value of cjj shown in the figure is u > j = 0.399 for 

k = 0.25, £ = 8 which is to be compared with o>j = 0.394 for k = 1.0, and £ = 0. 


Figure 6a shows only the curves for £ up to 6 * 10. This is because it is known, for the parameter 
values of Figure 6a, that the maximum for the pure symmetric modes occurs at £ = 8.8. This result is 

shown in Figure 6b. It was anticipated that the values of cjj will decrease further as £ increases beyond 

10 making additional curves unnecessary. As pointed out above, for Ri = 0.80 the maximum growth 
rates of the symmetric and conventional baroclinic modes are about the same. Stone [12] found for 
his inviscid analysis, that the maximum growth rates are equal for Ri = 0.95. The reduction in this 
value of Ri must be due to the inclusion of viscous and thermal effects in our model. Also all of the 
unstable modes were found to propagate with the mid-channel speed, U = 0 (i.e., for all cj; > 0, u> T = 
0). This result was found to hold for all the unstable modes presented in this study. 

Note that in Figure 6a the growth rate maximum does not occur for a pure symmetric wave; 
k at the maximum value is small but not zero. This oint will be elaborated upon further in the follow- 
ing discussion. To avoid confusion, symmetric instability is defined to include unstable modes for which 
k is small as well as zero. Further, for simplicity, conventional baroclinic instability (£ = 0) is referred to 
as zonal instability. In addition the term preferred mode will be used to mean the mode of maximum 
growth rate. 

Figure 7a for Ri = 0.50 shows the effect on ojj of a further decrease in Ri. The growth rates for 

the symmetric modes have continued to increase and are now considerably greater than those for the 
pure zonal modes. This figure shows only the curves for £ up to £ = 10. Again this is because, for the 
values of the parameters of this figure, the maximum of wj for the pure symmetric modes occurs at 

fi = 1 0.2. This result is shown in Figure 7b. For the data in Figure 7a the maximum growth rate < ccurs 
for k a 0.05 and £ = 10. Note again that the maximum does not occur for k = 0. This is clearly shown 
in Figure 7c which is a replot of the growth rate for the £ = 10 wave with a stretched coordinate. This 
means that the preferred mode of symmetric baroclinic instability exhibits a slight angle of inclination 
with the direction of the basic state flow. This angle is given approximately by k/£. An analogous result 
was found earlier by Busse and Chen [16] (see Section I) for the waves at maiginal stability. 

Figure 8a for Ri = 0.286 shows the effect on coj of decreasing Ri still further. This reduction 
leads to an increase in Wj for both k and £, but the overall functional dependence of on k and £ has 

not changed much from the results for Ri = 0.50. Again the preferred wave shows a slight inclination 
with respect to the basic state flow (k = 0.07, £ = 10) as is shown in Figure 8b. The magnitude of the 
inclination is a'most the same as that of Figure 7c. 

The results for another value of o are now presented. Figures 9 through 1 1 show «j versus k for 

three values of Ri, a = 5 and AT = -4.5. Figure 9, for Ri = 1.0, shows results that are qualitatively 
similar to those for o = 2, Ri = 0.93 (Fig. 5). Again pure symmetric modes possess some growth rate, 
but the dominant growth rate is found associated with the zonal modes. Figure 10a, for Ri = 0.50, shows 
results similar to those in Figure 7a. For this value of Ri, the symmetric modes are the preferred modes 
and the zonal modes have the lowest growth rates. Figure 10b, which is a stretched plot of the £ = 10 
curve in Figure 10a, shows that again the preferred mode is not purely symmetric but occurs for k as 
0.10. Figure 11a, for Ri = 0.135, shows results qualitatively similar to those in Figure 8a. These results 
are also similar to Figure 10a hut show a discernible increase in the growth rate values. Figure lib, 
which is a stretched plot of the £ = 10 curve in Figure 11a, shows that again the preferred mode is not 
purely symmetric but occurs for k as 0.10. Stone [12] (see Section I) found, for his inviscid analysis, 
that the preferred mode of instability changed at Ri * 0.25 from symmetric instability to Kelvin- 
Helmholtz instability. The results in Figure 1 la, for Ri < 0.25, show that symmetric waves are preferred. 
However, Stone’s conclusion was based on growth rate behavior for asymptotically large k; and, since 




growth rates for large k for this model were not computed, nothing definite can be said about such a 
transition for this analysis. It is possible, since viscous effects have been included, that this transition 
may occur only for very small E or not at all. Note that the earlier investigation of this model for pure 
symmetric modes (AF) revealed that the growth rates continue to increase as Ri is decreased through 
Ri < 0. 

Results for a < 1 are now presented. Figure 12a shows coj versus k for Ri = 0.50, a = 0.20 and 

AT = 0.50. Clearly, Ri = 0.50 is below the threshold value at which symmetric modes become the pre- 
ferred modes. A comparison of Figure 12a with Figures 7a and 10a shows that symmetric baroclinic 
instability is dominant at Ri = 0.50 for both a $ 1. A difference, however, in the sign of the angle of 
inclination of the symmetric waves with respect to the basic state flow occurs depending on whether 
o $ 1. Figure 12 b, which is again a replot of the £ = 8 curve in Figure 12a with a stretched coordinate, 
shows that the growth rate maximum occurs at k -0.11. A negative vah.e of k implies a change in 
direction of the x-component of the wavenumber vector. Thus, the angle of inclination of the symmetric 
waves with respect to the basic state flow has changed sign. This result, for the modes of maximum 
growth rate, is analogous to that found by Busse and Chen [ 1 6] for marginally stable waves. 

Attention is now turned to the results for a = 1, which were totally unexpected. Figures 13 
through 15 show cjj versus k for three values of Ri, a = 1.0 and AT = -0.20. Figure 13 for Ri = 0.80 

shows that the pure zonal waves are the preferred modes but that the symmetric modes also have sig- 
nificant growth rates. These results are qualitatively similar to the results for different values of a and 
similar values of Ri (Figs. 5 and 9). However, as Ri is reduced, the opposite general trend is found. 
Figure 14 for Ri = 0.50 shows that the growth rates of the symmetric modes have decreased significantly 
with respect to the zonal modes. Figure 15 for Ri = 0.30 shows that these relative growth rates have 
decreased still further. In fact they have vanished altogether. Busse and Chen [16] predicted, for mar- 
ginally stable waves for small k and a = 1, that the first-order correction to the value of Ri c is zero, 

indicating that the angle of inclination of the waves with respect to the basic state flow is also zero, 
rhis result does not agree with that found for the unstable waves for small k in this study. Examination 
of Figures 13 and 14 shows that the symmetric growth rate curves approach the ordinate with negative 
gradients implying a nonzero angle of inclination with the same sign as that for a < 1 results. 

Bus:e and Chen [16] showed that when zonal structure is allowed, Ri c increases above the pure 

symmetric value. Their analysis was limited to small values of k, but they pointed out that owing to 
this effect the range of Ri for which the instability occurs may be increased significantly beyond the 
limits observed by McIntyre [17], Since the numerical procedure is not limited to small k, it was decided 
to investigate this suggestion by extending the work of Busse and Chen to larger values of k. Figures 

16a and 16b show «j/£ as a function of Ri for several values of k for a = 2, AT = -1.5, e = 10'^, and 
for £ = 4 and 6, respectively. The choice of £ = 4 and 6 is based on previous work (AF) which showed 
maximum values of Ri c close to these values of £. In Figures 16 the value of Ri c for a specific value 

of k is the value of Ri at which the growth rate curve intercepts the Ri axis. Figure 1 7 is based on the 
results in Figure 16 and shows Ri c versus k for the two values of £. For both values of £ the qualitative 

functional behavior of Ri c with k is similar. At first Ri c increases with increasing k and rises above its 
value for k = 0. However, for still relatively small values of k, Ri c reaches a maximum vai te which is 
only about 10 percent greater than its value at k = 0, and then Ri c declines monotonically as k increases 

further. Thus, it has been shown, for the parameters of this computation, that the range of symmetric 
baroclinic instability is increased but not significantly when zonal structure is permitted. 
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IV. CONCLUSIONS 




t ^ - 






This work has been concerned with the three-dimensional baroclinic instability problem for Ri ~ 
0(1). The inviscid analyses of Stone [11,12,13] has been extended by adding viscous and thermal 
diffusion effects to the basic state and to the stability analysis and by satisfying realistic boundary condi- 
tions. Some pointed to Ri as the significant parameter and calculated the values of Ri at which the transi- 
tion between the conventional baroclinic instability, symmetric baroclinic instability, and the Kelvin- 
Heimholtz instability occur (see Section 1). Because of computing time limitations, it was not possible 
to perform a complete survey of these transitions for the model, but for 0 < Ri < 1.1, 0.2 < a < 5, 

E = 10'^, several flows were thoroughly studied. Although the results are not in quantitative agreement 
with Stone’s and several significant differences exist, the same general trend of the dominance of 
symmetric instability over conventional baroclinic instability as Ri decreases was found, except for a = 1. 
For a § 1 zonal instability dominates for Ri > 0.8, and symmetric instability dominates for Ri < 0.8. 
For a - 1 in the range 0.3 < Ri < 0.8, symmetric instability never dominates. 

It was found for a $ 1 that when the symmetric instability dominates, the preferred mode is not 
pure symmetric but has weak zonal structure. This means that the wave fronts are inclined at a small 
angle to the basic state eastward flow (x-direction) and that the wavenumber vector is inclined at the 
same angle to the north (y-direction). For a > 1, the vector has a component towards the east, and for 
a < 1 a component towards the west. Stone found the preferred mode to be truly symmetric. 

Busse and Chen [16] found, in the limit of small E and for small k, that the mode of symmetric 
instability at maximum Ri c has an angle of inclination with the basic state flow direction. This result is 

analogous to what was found for the symmetric modes of maximum growth rate. On the basis of their 
result, Busse and Chen speculated that the range of Ri for which symmetric instability occurs may be 
increased substantially over the range for k = 0. The results of this report, which are valid over a wide 
range of k, show that this is not so; only a small increase in Ri c occurs at small k. These results show 

that for values of Ri close to the values of Ri c for symmetric instability, the zonal instability has much 

larger growth rates and hence will dominate. This result reduces the significance of previous discussions 
on the behavior of Ri c . 

It has been shown that the nearly symmetric modes of maximum growth rate differ very slightly 
from the pure symmetric modes. For this reason it is unlikely that they represent a new type of 
instability. 


It should be recaiied that the conclusions of this report are based on calculations for which E 

3 

was kept at a constant value of 10 . Some modifications to these conclusions are anticipated if E is 

increased above 10 . However, based on our experience in AF, we are confident of no significant 
changes in the conclusions as E -*■ 0. 


It is of value to ask what the new results obtained in this study tell about the possibility of 
realizing symmetric baroclinic instability in the laboratory. The model is physically realistic because the 
full set of governing equations, including viscous and thermal diffusion, was used and realistic 
boundary conditions on the horizontal surfaces were satisfied. The model departs from laboratory reality 
in that a linear stability analysis rather than a nonlinear analysis was performed. Further, a horizontally- 
infinite layer is, of course, unrealizable, but a shallow and wide cylindrical annulus of fluid can be 
realized. Although the results of Section III indicate that the preferred mode of symmetric instability 
occurs with a large zonal wavelength and is not purely symmetric, this may not be so in finite cylindrical 
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geometry when the azimuthal wavelength exceeds the circumference. An apparatus has been constructed 
in which a shallow layer of liquid is held between two horizontally-mounted discs. Radial temperature 
gradients with different imposed vertical temperature differences can be maintained on the discs. The 
discs are made of sapphire material, so they are transparent allowing observation of the liquid. Sapphire 
is also a good thermal conductor allowing for an accurate temperature boundary condition. The total 
apparatus is rotated on a turntable. We have started a systematic search with 0 < Ri < 1 and are looking 
for convective rolls with radial sti-rture and weak or no azimuthal structure. 


1 

I 


♦ 


i 


11 



«• iAr • 



REFERENCES 




Pedlosky, J.: Geophysical Fluid Dynamics. Springer-Verlag, 1979. 

Chamey, J. G.: The Dynamics of Long Waves in a Baroclinic Westerly Current. J. Meteor., Vol. 
4, 1947, pp. 135-162. 

Eady, E. T.: Long Waves and Cyclone Waves. Tellus, Vol. 1, 1949, pp. 35-52. 

Lorenz, E. N.: The Nature and Theory of the General Circulation of the Atmosphere. World 

Meteorological Organization, 1967. 

Hide, R. and Mason, P. J.: Sloping Convection in a Rotating Fluid. Advances in Physics, Vol. 

24, 1975, pp. 47-100. 

Stone, P. H.: An Application of Baroclinic Stability Theory to the Dynamics of the Jovian 

Atmosphere. J. Atmos. Sci., Vol. 24, 1967, pp. 642-652. 

Bennets, D. A. and Hoskins, B. J.: Conditional Symmetric Instability a Possible Explanation 

for Frontal Rainbands. Quart. J. Roy. Met. Soc., Vol. 105, 1979, pp. 945-962. 

Emanuel, K. A.: Inertial Instability and Mesoscale Convective Systems. Part I: Linear Theory of 
Inertial Instability in Rotating Viscous Fluids. J. Atmos. Sci., Vol. 36, 1979, pp. 2425-2449. 

Solbeig, H.: Le Mouve.nent d’ inertie de i’atmosphere stable et son role dans le theorie des 

cyclones. Union Geodesique et Geophysique Internationale Vlieme Assembiee. Edinburgh, Vol. 
II, 1936, pp. 66-82. 

Kuo, H. L.: Forced and Free Axially Symmetric Convection Produced by Differential Heating 
in a Rotating Fluid. J. Meteor., Vol. 13, 1956, pp. 521-527. 

Stone, P. H.: On Ncn-Geostrophic Baroclinic Stability. J. Atmos. Sci., Vol. 23, 1966, pp. ^90- 
400. 

Stone, P. H.: On Non-Geostrophic Baroclinic Stability: Part II. J. Atmos. Sci., Vol. 27, 1970, 
pp. 721-726. 

Stone, P. H.: Baroclinic Stability Under Non-Hydrostatic Conditions. J. Fluid Mech., Vol. 45, 
1971, pp. 659-671. 

Stone, P. H., Hess, S„ Hadlock, R., and Ray, P.: Preliminary Results of Experiments with 
Symmetric Baroclinic Instability. J. Atmos. Sci., Vol. 26, 1969, pp. 991-996. 

Hadlock, R. K., Na, J. Y., and Stone, P. H.: Direct Thermal Verification of Symmetric Baro- 
clinic Instability. J. Atmos. Sci., Vol. 29, 1972, pp. 1391-1393. 


Busse, F. H. and Chen, W. L.: On the (Nearly) Symmetric Instability. J. Atmos. Sci., Vol. 38, 
1981, pp. 877-880. 



McIntyre, M. E.: Diffusive Destabilization of the Baroclinic Circular Vortex. Geophys. Fluid 

Dyn., Vol. 1, 1970, pp. 19-57. 

Antar, B. N. and Fowlis, W. W.: Symmetric Baroclinic Instability of a Hadley Cell. J. Atmos. 
Sci., VoL 39, 1982, pp. 1280-1289. 

Antar, B. N. and Fowlis, W. W.: Baroclinic Instability of a Rotating Hadley Cell. J. Atmos. 

Sci., Vol 38, 1981, pp. 2130-2141. 


Conte, S. D.: The Numerical Solution of Linear Boundary Value Problems. S.I.A.M. Rev., 

Vol. 8, 1966, pp. 309-320. 



Ri c 1.0 


o 


o 


ol 1 1 

0.1 1.0 10 


Figure la. The critical Richardson number as a function of the Prandtl number, according to 
McIntyre [17] and according to Antar and Fowlis [18]. 



Figure 1 b. Marginal stability curves as functions of E and Ri for three values of o and 
for £ = 4. The dashed curves are McIntyre’s results. 
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Figure 3a. The basic state zonal velocity profile as a function of height according 
to equation (6) for three values of the Ekman number. 
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Figure 4. The temperature profile for E = 10'^ and three values of the Prandtl number; 
o=l, AT = -0.2, and e = 1.6; o = 2, AT = -1.5, and e = 4.0; 
o = 5, AT = -4.5, and e = 4.0. 
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Figure 5. The growth rates as a function of k for Ri = 0.93, o = 2, AT = -1.5, E = 10'' 

and selected values of £. 
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Figure 1 la. The growth rates as a function of k for Ri = 0.135, a - 5, 
AT = -4.5, E = 10* 3 and selected values of £. 



Figure 1 lb. The growth rates as a functic A k for £ = 10, Ri = 0.135, a = 5, 

AT = -4.5, and E = 10' 3 . 
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Figure 16a. 
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growth rates as a function of k for Ri = 0.30, o = 1.0, AT = -0.2 
E = 10* 3 and selected values of &. 



: normalized growth rates as a function of Ri for £ = 4, a = 2, 
AT = -4.5, E = 10‘ 3 and selected values of k. 
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Figure 1 7. The critical Richardson number, Ri c as a function of k for fi = 4 and 6, 

a = 2, AT = -1.5, and E = 10' 3 . 
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